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ABSTRACT 

We have extended our previous study of the lattice QCD spectrum with 2 flavors of staggered 
dynamical quarks at 6/(/^ = 5.6 and aniq = 0.025 and 0.01 to larger lattices, with better 
statistics and with additional sources for the propagators. The additional sources allowed us 
to estimate the A mass and to measure the masses of all mesons whose operators are local 
in time. These mesons show good evidence for flavor symmetry restoration, except for the 
masses of the Goldstone and non-Goldstone pions. PCAC is observed in that oc m^, and 
is estimated. Use of undoubled lattices removes problems with the pion propagator found 
in our earlier work. Previously we found a large change in the nucleon mass at a quark mass 
of aniq = 0.01 when we increased the spatial size from 12 to 16. No such effect is observed 
at the larger quark mass, aniq = 0.025. Two kinds of wall source were used, and we have 
found difficulties in getting consistent results for the nucleon mass between the two sources. 



1 INTRODUCTION 



Calculations of liadron spectroscopy remain an important part of nonperturbative studies of 
QCD using lattice methods. (For reviews of recent progress in this field, see Ref. [1].) We 
have been engaged in an extended program of calculation of the masses and other parameters 
of the light hadrons in simulations that include the effects of two fiavors of light dynamical 
quarks. These quarks are realized on the lattice as staggered fermions. We have carried 
out spectrum calculations with lattice valence quarks in both the staggered and Wilson 
formulations. 

These simulations are performed on f6"^x32 lattices at lattice coupling /3 = 6/(/^ = 5.6 
with two masses of dynamical staggered fermions, aniq = 0.025 and aniq = 0.01. These are 
the same parameter values as we used in our first round of simulations [2]. However, the 
first set of simulations had two known inadequacies. The first was that most of our runs 
were carried out on lattices of spatial size 12"^. A short run on 16'* lattices with dynamical 
quark mass 0.01 showed that the 12"* lattices were too small: baryon masses fell by about 
fifteen per cent on the larger lattice compared to the smaller one. Thus, it was important 
to investigate finite size effects for aniq = 0.025. We also felt the need for more statistics on 
the aniq = 0.01 system for lattices of spatial size 16"^. 

Second, nearly all of our earlier running was done on lattices of size 12"* or 16"*; these 
lattices were doubled (or quadrupled) in the temporal direction to 12"^ X 24 (or 12"^ X 48) or 
16"^ X 32 for spectroscopy studies. Doubling the lattice introduced structure in the propagators 
of some of the particles: the pion effective mass, in particular, showed peculiar oscillatory 
behavior as a function of position on the lattice. This behavior was almost certainly due 
to doubling the lattice [3] and the best way to avoid this problem is to begin with a larger 
lattice in the temporal direction. Because of these difficulties, mass estimates from such 
doubled lattices are suspect. This is seen when comparing the masses obtained from the 
doubled or quadrupled 12"* lattice with those from the 12"^ X 24 lattice in our previous work. 
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In our work on smaller lattices, only one kind of source was used, the so-called "corner" 
source. In these simulations we include a second kind of source (in fact 3 sources) which 
enables us to measure the A mass, as well as the nucleon. Furthermore, with these new 
sources, we are able to measure masses of all mesons created by operators which are local 
in time, and correspond to strictly local continuum operators (local quark bilinears with no 
derivatives). This allows us to study the extent to which flavor symmetry, which is broken 
by the staggered lattice, is realized at this lattice spacing. 

Some of the results described here have been presented in preliminary form in Ref. [4]. 
Studies with Wilson valence quarks which complement the results presented here have been 
published [5] as have studies of Coulomb gauge wave functions [6]. In addition, we are 
preparing a paper on glueballs and topology. In Sec. 2 we describe our simulations and in 
Sec. 3 we give our results and conclusions. 

2 THE SIMULATIONS 

Our simulations were performed on the Connection Machine CM-2 located at the Supercom- 
puter Computations Research Institute at Florida State University. 

We carried out simulations with two flavors of dynamical staggered quarks using the 
Hybrid Molecular Dynamics algorithm [7]. The lattice size was 16"^ X 32 sites and the lattice 
coupling /3 = 5.6. The dynamical quark masses were aniq = O.OI and 0.025. A total of 
2000 simulation time units (with the normalization of Ref. [2]) was generated at each value 
of the quark mass, after thermalization. The aniq = O.OI run started from an equilibrated 
16'* lattice of our previous runs on the FTA-IO, which was doubled in the time direction 
and then re-equilibrated for 150 trajectories. The aniq = 0.025 run was started from the 
last configuration of the smaller mass run, and then thermalized for 300 trajectories. For 
aniq = O.OI, we recorded lattices for the reconstruction of spectroscopy every 5 HMD time 
units, for a total of 400 lattices. At aniq = 0.025, lattices were stored every 10 time units 
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for a total of 200 lattices. 

For our spectrum calculation, we used periodic boundary conditions in the three spa- 
tial directions, and antiperiodic boundary conditions in the temporal direction. To calculate 
hadron propagators, we fixed the gauge in each configuration to the lattice Coulomb gauge 
using an overrelaxation algorithm[8], and used sources for the quark Green functions which 
spread out in space uniformly over the spatial simulation volume and were restricted to a 
single time slice ("wall" sources [9]). Our inversion technique was the conjugate gradient 
algorithm, using a fast matrix inverter written in CMIS (a low level assembler for the CM-2) 
by C. Liu[10, 11]. 

In this work we used two kinds of wall sources. The first of these consisted of a 1 
in a selected color component at each site of the source time slice where the x, j/, and z 
coordinates were all odd. In other words, the source was restricted to a single corner of each 
2^ hypercube. This is the same source as used in our previous work, and we will refer to it 
as the "corner" source or C. 

In addition to this corner source, we also used a triplet of wall sources. Following 
Gupta et al. [12] we defined an "even" source which takes the value +1 on every site of the 
source time slice, and an "odd" source which is +1 on all the even (space odd) sites on the 
source time slice and -1 on all the odd (space even) sites on this time slice. For definiteness, 
in this paragraph we take the source time slice to be t = 1. These sources allowed calculation 
of the A propagator and propagators for some of the local and non-local mesons. (We use 
the A propagator corresponding to a point sink where the three quarks are displaced by one 
unit in the x, j/, and z directions, respectively, from the origin of the unit cube [13].) The 
third source we used was what we call a "vector" source. This source is -|-1 on all sites on 
the source time slice that have an even y coordinate and -1 for those whose y coordinate is 
odd. With these three sources we are able to calculate meson propagators for all 20 meson 
representations of the time slice group which are local in time [14]. In addition, we have 
calculated the propagator for a local nucleon, and the A discussed above, from the "even" 
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source quark Green functions. We will refer to this triplet of sources as EOV. 

For the mesons, we averaged propagators computed from six sets of wall sources at 
time slices 1, 2, 3, 17, 18 and 19. Propagators from three consecutive time slices were needed 
for a separate study of glueball to qq correlations. For the baryons, we used four wall sources 
at time slices 1,9, 17 and 25. 

Finally, for comparison, we also measured the hadron propagators from a point source. 
This calculation was performed "on line" every time unit, for a total of 2000 measurements 
for each quark mass. 

3 RESULTS 

3.1 Doubling effects on the pion propagator 

In our previous, work we used 12"* and IG"* lattices doubled or quadrupled in the time direction 
for computing propagators. We found irregularities in the effective mass as a function of 
distance from the source for the pion. (The effective mass is the mass obtained by fitting 
with zero degrees of freedom to points in the propagator centered at some distance.) For 
the pion propagator, where we fit to a simple exponential plus the piece from periodicity, 
the effective mass at distance T is obtained from the two points in the propagator at time 
distances T — | and T + |. For other particles, where we use four parameter fits, with one 
particle of each parity, the effective masses are obtained from four successive distances in 
the propagator ( T— |, |, r+|, r+|). Since the location of these features seemed 
to depend on the lattice size before doubling, we tentatively ascribed them to effects of the 
doubling [2]. A simple analytic model of a doubled lattice showed similar features [3]. In the 
current work, we generated configurations on a 16"^ X 32 lattice, and did not double in the 
time direction when computing the hadron spectrum. The pion propagator is much better 
behaved. We show the new results for the pion effective mass together with our previous 
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Figure 1: Pion effective mass versus distance for aniq = 0.025. 

results in Figs. 1 and 2, for aniq = 0.025 and 0.01 respectively. Note that one of the sources 
used in the present work, the "corner" wall source, is identical to the source used in the 
previous work. We see that the pion effective mass in the current work is the same for the 
two sources, and is much smoother than on the doubled lattices. 

In Figs. 3 and 4, we show the effective mass plots for the p. Again, we notice that 
they are relatively flat, in contrast with the work using doubled lattices. [2] 
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Figure 2: Pion effective mass versus distance for aniq = 0.01. 

3.2 Best estimates for masses 

Hadron masses were estimated by making correiated fits to the average propagator[15]. To 
reduce the effects of autocorreiations in simuiation time, propagators on severai successive 
iattices were averaged together before computing the covariance matrix. For exampie, we 
most commoniy biocked 8 iattices together for aniq = 0.01, or 40 time units, since we 
measured every five time units. For aniq = 0.025 we typically blocked together 4 lattices or 
again 40 time units. 

To display the fits we use figures in which the symbol size is proportional to the 
confidence level of the fits. The symbol size in the keys corresponds to a confidence level of 
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Figure 3: Rho effective mass versus distance for aniq = 0.025. 

0.5. We plot the fits as a function of the minimum distance used in the fit. To show how the 
fit quahty varies with distance from the source, we plot the fits with two degrees of freedom. 
Such fits for the p masses are displayed in Figs. 5 and 6 and for the nucleon in Figs. 7 and 8. 

Tables 1-4 give our estimates for the hadron masses. In the continuum, all 15 compo- 
nents of the TT multiplet should be degenerate, as should all 15 components of the p multiplet. 
(Although we have only two fiavors of quarks in internal lines, the external quark lines have 
four quark fiavors. Hence, in the continuum limit hadrons form multiplets of fiavor 5^7(4).) 
When using staggered quarks on the lattice, fiavor symmetry is explicitly broken, and each 
continuum fiavor multiplet is broken down into irreducible representations of the discrete 
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Figure 4: Rho effective mass versus distance for aniq = 0.01. 

symmetries of the lattice action restricted to a given time slice [14]. Full flavor symmetry 
should be restored in the continuum limit. The extent to which this symmetry is restored 
at a flnite lattice spacing gives us some indication as to whether our lattice spacing is small 
enough (/3 large enough) to adequately approximate the continuum limit. In Figs. 9 and 10, 
we plot the tt and p masses, respectively, from Tables 1-4 for the different representations 
of the time slice group accessible using the FOV sources. We notice that flavor symmetry 
appears to be good to a few percent for the p multiplet. For the tt sector there is approxi- 
mate degeneracy for the non-Goldstone pions, but the mass of the Goldstone pion still lies 
signiflcantly below that of the rest ((m~ — m^)/m^ 0.3 for aniq = 0.01. This should not 
surprise us since, for the quenched approximation, deflnitive evidence for the restoration of 
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Figure 5: Fits to the p mass for aniq = 0.025. The size of the points is proportional to the 
confidence level of the fits. 

the mass degeneracy between the Goldstone and non-Goldstone pions has only been claimed 
for /3 = 6.5 [16]. In Sec. 3.5, we shall indicate that our /3 (5.6) is more comparable with a 
quenched system at /3 5.95. The squared Goldstone pion mass is very nearly proportional 
to niq (see Sec. 3.4). The other pion masses do not extrapolate to zero with niq. For ex- 
ample, at aniq = 0.025 the mass ratio m~/m^ = 1.223(5), while at aniq = 0.01 this ratio is 
1.306(11). This contrasts with a four fiavor study by the MTc collaboration[17], in which 
m~ appears to be proportional to niq. 

Fig. 11 gives the "Fdinburgh" plot of mj^/nip against mj^/nip for the results of Ta- 
bles 1-4. Fig. 12 is the "Boulder" plot for the — A mass splitting. Both plots are roughly 
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Figure 6: Fits to the p mass for aniq = 0.01. The size of the points is proportional to the 
confidence level of the fits. 

what one would expect for these values of the quark mass. 

3.3 Finite size and source effects on the hadron masses 

In our earlier work, we found a large change in the nucleon mass with quark mass aniq = 0.01 
when the spatial lattice size was increased from 12 to 16. With our new results, we can 
examine this in more detail, as well as extend the study to aniq = 0.025. 

For aniq = 0.01, the nucleon mass on a 12"^ X (12 X 2) lattice was estimated to be 
0.848(11) while that on a 12"^ X 24 was found to be 0.815(13), the difference probably being 
an effect of the doubling. On the 16"^ X (16 X 2) this had fallen to 0.770(8). In the data of 
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Figure 7: Fits to the nucleoli mass for aniq = 0.025. The size of the points is proportional 
to the confidence level of the fits. 



Table 1 for the 16"^ X 32 lattice we find that for the "corner" source (which is identical to the 
source used on the smaller lattices) the value is 0.748(4) again lower than the doubled case. 
Thus, we have further evidence for the finite volume effect reported in [2] and also seen by 
[18] for the nucleon mass. At aniq = 0.025, the nucleon mass on a 12"^ X (12 X 2) lattice was 
0.982(9), while that for a 16"^ X 32 lattice (Table 3) is 0.981(8). Thus, it would appear that 
even a 12"^ box is adequate to hold a nucleon at aniq = 0.025 with no appreciable finite size 
effects. 

For the mesons, we find good agreement between the masses on 12"^ X 24, 16"^ X (16 X 2) 
and the new results on 16"^ X 32 lattices, for both quark masses. The 12"^ X (12 x2) p effective 
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Figure 8: Fits to the nucleoli mass for aniq = 0.01. The size of the points is proportional to 
the confidence level of the fits. 

masses showed no clear plateau, and this is refiected in the other fits, so that the mass 
values were unreliable. The observed undulating behavior of the tt effective masses on the 
12^ X (12 X 2) lattice refiects itself in the more general fit. (The 16^ X (16 X 2) lattice 
shows similar problems.) Within these ambiguities, the new results are in good agreement 
with those for smaller lattices. Hence we may conclude that there are no significant finite 
size effects in the meson masses for spatial boxes with volumes > 12"^ for quark masses 
aniq > 0.01 at /3 = 5.6. 

Now let us discuss the effects of the two different types of source. For mesons, there 
is only one wall source, and one point sink corresponding to each component of each irre- 
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Figure 11: Edinburgh plot. The diamonds are our resuhs from the previous simulations, 
and the squares are the new results. Error bars are the statistical errors only, and do not 
include uncertainty based on the choice of source or fitting range. (We use the even-odd 
source results here.) 

ducible representation of the time slice group [14]. This means that for those representations 
occurring in both the "corner" and EOV wall sources we can expect to get the same results 
in both cases. This is well born out by the masses of Tables 1-4. 

Eor the nucleon, we use a local sink which projects the 8 representation of the time 
slice group. The "corner" source produces only one baryon representation, the local 8 
representation. The "even" source, on the other hand produces all baryon representations, 
and in particular 5 copies of the 8 representation. Only one of these 8's is local; the other 
4 have quarks on more than 1 vertex of the unit cube. The local point sink has overlap 



14 



0.15 



0.10 



0.05 



0.00 



T 1 1 \ ^ \ \ 1 r 

o This work 
X HEMCGC 1990 
□ APE |8 = 5.7 
o APE |8 = 6.0 




0.0 0.1 0.2 0.3 

(m -mj/(3m +m^) 



Figure 12: Comparison of baryon and meson hyperfine splitting. The two circles show the 
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the line interpolating between them is a simple quark model. The APE data is from a 
quenched simulation [19]. 



with to all 5 of these octets, each of which will, in general, have different couplings to the 
allowed baryon states. For this reason the nucleon propagator for the "even" source can 
be rather different from that for the corner source. That this is so is illustrated by looking 
at the effective mass plots (Figs. 13 and 14) for the 2 different nucleon propagators. At 
aniq = 0.025, the effective masses for the "corner" source lie consistently higher than those 
for the "even" source. Since it is difficult to find strong evidence for a plateau in this data 
(at least not for the "corner source") the problem could well be that the plateau starts just 
as the signal/noise ratio starts to worsen. In any case, our best fits (Table 3) are within 
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Figure 13: Nucleoli effective mass versus distance for aniq = 0.025. 



2 standard deviations of one another, and can thus be considered to be consistent. For 
aniq = 0.01, the effective masses for the "corner" nucleon again lie consistently above those 
for the "even" source. However, in the graph of Fig. 14, one notes that the effective masses 
for the 2 sources appear to be coming together for T > 7.5. If this is correct, the reason for 
the discrepancy between the two estimates of the nucleon mass is that our fitting criterion 
favors the false plateau 3.5 < T < 6.5 in the nucleon effective mass plot. If we had better 
statistics we would presumably find the true plateau. 

Finally let us comment on the point source fits as compared with the wall fits. Only 
in the case of the tt do these point source fits have the quality of the wall fits. The tt masses 
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Figure 14: Nucleoli effective mass versus distance for aniq = 0.01. 

obtained from the point and wall sources are in excellent agreement. For the other particles, 
the rapid decrease of the point source propagators with increasing T due to contamination 
with higher mass excitations produces mass estimates that tend to be high and at the very 
least have much larger errors than the wall results. The rapid decrease of effective masses 
with T for the point sources makes the evidence that these reach a plateau before the signal 
is lost less compelling than in the case of wall sources. Their main virtue is that their mass 
predictions give an upper bound on the particle mass. This is no great advantage if the 
bound is too large or has a very large uncertainty (error). 

Fventually, lattice QCD should provide masses for excited state hadrons as well as 
for the lowest mass particle for given quantum numbers. Identification of excited states is 
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Figure 15: Goldstone pion fits including excited state masses. The symbol size is propor- 
tional to the confidence level, with the symbol size used in the legend corresponding to 50% 
confidence. We show results for aniq = 0.025 (a) and 0.01 (b). The octagons, squares and 
bursts correspond to two particle fits, both particles having the same parity, with corner, 
EOV and point sources respectively, while the diamonds are from one particle fits to the 
corner wall source. 

probably easiest in the pion channnel, since the small ground state mass means that the 
excited state is probably well separated in mass from the ground state. In Fig. 15, we show 
fits to the pion from Fuclidean time range Tj^i^ to 16, including both the ground state and 
excited state masses for the two particle fits (both with the same parity). In these graphs 
the symbol size is proportional to the confidence level of the fits. At aniq = 0.025 we see 
consistent results for the excited state mass, independent of Tmin up to the point where the 
excited state is no longer needed in the fit, and independent of whether we use the FOV or 
the corner wall source. For the point source, we see that including two particles is still not 
sufficient to get good fits with small Tj^i^. Unfortunately, for aniq = 0.01 the results are not 
nearly as good, although this is one of the cases for which we have no good fits even with 
large T^m- 
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Figure 16: Squared pion masses versus quark mass. The octagons are for the Goldstone 
pion, and the squares for the other pointhke pion, the vr. The dashed and dotted hues are 
extrapolations to zero quark mass, and the horizontal lines on the left side are one standard 
deviation limits on this extrapolation. (In reality, we do not expect the vr mass squared to 
be vanish linearly with niq for small niq.) 



3.4 PCAC 



PCAC predicts that oc niq. In Fig. 16, we plot against niq. For definiteness we have 
chosen the 4 parameter FOV estimate for the pion mass in each case. We obtain 



0.0013(9) + 6.98(5)m„. 
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The intercept is only 1.4(7 from zero. Thus, this simple PC AC relationship appears to be 
well satisfied. We therefore can make use of the the more precise relationship 

for niq sufficiently small, to extract an estimate for To finess the question of perturbative 
subtractions for (V'V') which are known to remove most of its apparent mass dependence, 
we linearly extrapolate it to niq = where no such subtraction is necessary. Our measured 
values for a^{il)il)) were 0.11223(46) at am, = 0.01 and 0.21398(34) for am, = 0.025. This 
gives a^{t/jt/j) = 0.04440(80) at m, = 0. Since (V'V') measured with 4 fermion fiavors, we 
must multiply it by [Nj = 2)/4 = 0.5 before inserting it into the above equation giving 

af^ = 0.0564(5) . 

Estimating a by setting the p mass linearly extrapolated to m, = to its experimental value 
we find a = l/(1.80(2)GeV) (from the EOV p). This gives us 

A = 102(2)Mey 

as compared with the experimental value 93MeV . Note that the error we have quoted 
represents only the statistical error. Just taking into account the systematic uncertainty of 
choosing which p masses to use, indicates that the error estimate should almost certainly 
be at least twice what we have quoted. With this in mind, and remembering that we really 
have little justification for linearly extrapolating our p masses to m, = 0, we consider this 
value to be quite good. 

3.5 Comparison between quenched and full QCD 

In addition to simulating with 2 fiavors of dynamical quarks, we also estimated the hadron 
spectrum for quenched QCD on a 16"^ X 32 lattice. Here the aim was to compare the 
spectrum from quenched QCD with that of full QCD. For this reason, we chose /3 values 
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for the quenched runs which we beheved to be in close correspondence with /3 = 5.6 for full 
(2 flavor) QCD. Two values of /3 were chosen so that interpolation to a requisite /3 might 
be possible. We chose /3 = 5.85 and /3 = 5.95 as our two values. The hadron masses for 
these /3's at aniq = 0.01 and aniq = 0.025 are given in Tables 5-6. A cursory comparison of 
the masses in these tables with those for the full theory (Tables 1-4) shows why these two 
/3 values were chosen. At /3 = 5.85 the masses (in lattice units) of the Goldstone pions are 
very close to their values in the full theory for both quark masses. At /3 = 5.95, on the other 
hand, the masses of the p and nucleon at aniq = 0.01 are close to the values in full QCD. 

What is immediately clear from these results is that a simple shift in the coupling 
constant (i.e., in /3) is inadequate to reproduce the whole effect of including dynamical 
quarks as some have suggested. However, a /3 shift combined with a renormalization of 
the bare quark mass can bring the spectra into reasonable agreement. We find that the 
quenched spectrum at /3 = 5.95 is brought into reasonable agreement with the dynamical 
quark spectrum at /3 = 5.6 if we increase the bare masses in the quenched case by a factor 
of 1.16. The reason we must scale all quark masses by the same factor is because both the 
quenched and the dynamical Goldstone pions appear to obey PC AC. The new Goldstone 
pion masses are obtained using PCAC from those at aniq = 0.01 and aniq = 0.025. The non- 
Goldstone pion mass is obtained by noting that the difference between the non- Goldstone 
and Goldstone pion masses depends only weakly on the mass. Rho and nucleon masses 
for the "renormalized" masses are obtained from those at aniq = 0.01 and aniq = 0.025 by 
linear interpolation/extrapolation. The comparison between these quenched and unquenched 
masses is exhibited in Tables 7-8. These results might have been improved still further if we 
had varied /3 in the neighborhood of /3 = 5.95. The mass differences between the full QCD 
and quenched masses in Tables 7-8 are larger than can be attributed to statistics alone, but 
are probably consistent with the systematic errors due to choices of fits and differences in 
the finite size/lattice spacing errors between the two theories. It therefore remains to be 
seen whether there are significant differences in the infinite volume continuum theories. 
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arriq = 0.01 



particle 


source 


range 


mass 


error 


1 d.o. f. 


confidence 


parameters 


TT 


point 


12-16 


0.2681 


0.0010 


0.13 


0.95 


2 


TT 


EOV 


7-16 


0.2667 


0.0008 


1.90 


0.55 


2 


TT 


EOV 


1-16 


0.2673 


0008 


1.50 


0.12 


4 


TT 


C 


13-16 


0.2667 


0.0015 


1.45 


0.24 


2 


TT 


C 


4-16 


0.2700 


0.0012 


1.33 


0.20 


4 


TT 


point 


9-16 


0.3899 


0.0190 


1.27 


0.28 


4 


TT 


EOV 


6-16 


0.3500 


0.0026 


0.885 


0.52 


4 


TT 


C 


9-16 


3553 


0039 


0.20 


0.94 


4 




EOV 


5-16 


0.3474 


0.0014 


1.99 


0.30 


2 


^3(1) 


EOV 


6-16 


0.3694 


0.0032 


0.77 


0.61 


4 


V / 


EOV 


7-16 


0.3703 


0.0019 


1.25 


0.19 


2 


V / 


EOV 


8-16 


0.3842 


0.0034 


1.30 


0.26 


4 


7rf3) 


EOV 


7-16 


0.3831 


0.0021 


0.89 


0.52 


2 


77' / TT 


EOV 


6-16 


0.3952 


0.0036 


0.44 


0.88 


4 


r 


D 01 n t, 


9-16 


0.492 


0.038 


0.65 


0.63 


4 





EOV 


3-16 


0.5133 


0.0022 


1.11 


0.35 


4 





C 


8-16 


0.5085 


0.0050 


1.54 


0.17 


4 


n 
H 


D 01 n t, 


8-16 


0.476 


0.035 


0.31 


0.91 


4 



r 


EOV 


4-16 


0.5152 


0.0032 


1.74 


0.07 


4 


n 
H 


C 


10-16 


0.4918 


0.0091 


0.40 


0.74 


4 




EOV 


2-16 


5206 


0025 


1.31 


0.21 


4 




EOV 


2-16 


0.5184 


0.0030 


0.85 


0.59 


4 




EOV 


4-16 


0.5207 


0.0025 


0.64 


0.76 


4 




EOV 


3-16 


0.5205 


0.0024 


1.13 


0.34 


4 




EOV 


5-16 


0.5173 


0.0035 


1.05 


0.39 


4 


F V"/ 


EOV 


2-16 


0.5180 


0.0032 


0.92 


0.52 


4 


F \/^J 


EOV 


3-16 


0.5221 


0.0019 


0.88 


0.55 


4 


F \/^J 


EOV 


5-16 


0.5134 


0.0038 


1.29 


0.24 


4 




EOV 


2-16 


0.5229 


0.0022 


0.75 


0.69 


4 


r V / 


EOV 


3-16 


0.5186 


0.0031 


0.73 


0.70 


4 


N 


point 


9-14 


0.738 


0.086 


0.47 


0.63 


4 


N 


EOV 


7-15 


0.720 


0.006 


0.36 


0.88 


4 


N 


C 


10-15 


0.696 


0.027 


2.00 


0.37 


4 


N 


C 


0-15 


0.727 


0.008 


0.31 


0.96 


8 


N' 


point 


NP 


NP 


NP 


NP 


NP 


4 


N' 


point 


7-13 


1.209 


0.087 


1.21 


0.30 


4 


N' 


EOV 


6-15 


0.948 


0.066 


0.37 


0.90 


4 


N' 


EOV 


0-15 


0.948 


0.025 


0.31 


0.96 


8 


N' 


C 


3-15 


0.904 


0.009 


1.78 


0.06 


4 


A 


EOV 


4-15 


0.850 


0.008 


0.39 


0.93 


4 


A' 


EOV 


6-15 


1.031 


0.065 


0.47 


0.83 


4 



Table 1: Hadron masses aniq = 0.01. Notation: superscript is dimension of representation 
of time slice group; number of links in parenthesis; tilde (~) state has extra 70; notation 
abbreviated when unambiguous. "NP" indicates no plateau was found in the mass fits. 
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aniq = 0.01 


particle 


source 


range 


mass 


error 


X /d.o.f. 


coniiaence 


parameters 


TT 


bO V 


1-16 


0.893 


0.021 


1.50 


0.12 


4 


TT 


c 


4-16 


0.578 


0.063 


1.33 


0.20 


A 

4 


TT* 


point 


5-16 


0.789 


0.033 


1.87 


0.06 


4 


Jo/ao 


point 


10-16 


0.547 


0.015 


1.61 


0.19 


4 


Jo/ao 


bjU V 


b-lb 


(J. 514 


(J.UUo 


n on 


n IX 


A 

4 


Jo/ao 


c 


7-16 


0.505 


0.013 


0.78 


0.59 


A 

4 


3/1 \ 

«o(l) 


bO V 


6-16 


0.615 


0.014 


0.77 


0.61 


4 


3 /o\ 

«o(2) 


bO V 


6-16 


0.615 


0.019 


1.53 


0.15 


4 


ao(3) 




6-16 


0.645 


0.020 


0.44 


0.88 


4 


ai 


point 


8-16 


0.683 


0.097 


0.31 


0.91 


4 


ai 




5-16 


0.700 


0.011 


1.88 


0.06 


4 


ai 


C 


6-16 


0.744(.'^) 


0.020 


2.71 


0.008 


4 




rjU V 


o— ib 


U.byo 


U.UU ( 


U.yo 


U.oi 


4 




EOV 


6-16 


0.712 


0.013 


1.18 


0.31 


4 


a?(2) 


EOV 


5-16 


0.655 


0.018 


0.57 


0.80 


4 


«?(2) 


EOV 


3-16 


0.701 


0.004 


1.27 


0.24 


4 


ai(3) 


EOV 


4-16 


0.726 


0.011 


0.80 


0.62 


4 




point 


7-16 


0.818 


0.135 


0.95 


0.46 


4 


w 


EOV 


3-16 


0.686 


0.008 


1.11 


0.35 


4 


w 


C 


6-16 


0.775(?) 


0.042 


1.86 


0.08 


4 


m) 


EOV 


2-16 


0.719 


0.007 


1.31 


0.21 


4 




EOV 


5-16 


0.739 


0.015 


0.63 


0.76 


4 


/^i/fei 


EOV 


5-16 


0.732 


0.019 


1.05 


0.39 


4 




EOV 


2-16 


0.717 


0.004 


1.05 


0.40 


4 


W{3) 


EOV 


2-16 


0.711 


0.006 


0.75 


0.69 


4 



Table 2: Hadron masses aniq = 0.01. Notation: superscript is dimension of representation 
of time slice group; number of links in parenthesis; tilde (~) state has extra 70; notation 
abbreviated when unambiguous. The "?" denotes cases where none of the fits were good. 
The TT* is an excited state in the pion channel. 
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anig = 0.025 


particle 


source 


range 


mass 


error 


xVd.o.f. 


confidence 


parameters 


TT 


point 


13-16 


0.4188 


0.0005 


1.67 


0.17 


2 


TT 


point 


4-16 


0.4190 


0.0005 


1.68 


0.09 


4 


TT 


EOV 


10-16 


0.4185 


0.0009 


0.90 


0.48 


2 


TT 


EOV 


1-16 


0.4193 


0.0007 


1.13 


0.33 


4 


TT 


C 


9-16 


0.4185 


0.0006 


0.63 


0.70 


2 


TT 


C 


1-16 


0.4192 


0.0006 


0.45 


0.94 


4 


TT 


point 


9-16 


0.5120 


0.0061 


1.43 


0.22 


4 


TT 


EOV 


8-16 


0.5106 


0.0018 


1.05 


0.39 


4 


TT 


C 


10-16 


0.5089 


0.0027 


0.63 


0.59 


4 


7r3(l) 


EOV 


9-16 


0.5098 


0.0012 


0.72 


0.64 


2 


^3(1) 


EOV 


8-16 


0.5347 


0.0025 


1.20 


0.30 


4 


7r3(2) 


EOV 


9-16 


0.5297 


0.0017 


0.42 


0.87 


2 


^3(2) 


EOV 


9-16 


0.5451 


0.0028 


0.96 


0.43 


4 


7r(3) 


EOV 


9-16 


0.5431 


0.0018 


0.69 


0.66 


2 


r]'/Tr 


EOV 


6-16 


0.5507 


0.0020 


0.53 


0.82 


4 


P 


point 


10-16 


0.6243 


0.0127 


0.40 


0.75 


4 


P 


EOV 


7-16 


0.6396 


0.0028 


1.33 


0.24 


4 


P 


C 


8-16 


0.6396 


0.0037 


1.00 


0.42 


4 


P 


point 


NP 


NP 


NP 


NP 


NP 


4 


P 


EOV 


8-16 


0.6471 


0.0043 


0.77 


0.57 


4 


P 


C 


6-16 


0.6437 


0.0037 


0.60 


0.76 


4 


(jj/p 


EOV 


8-16 


0.6507 


0.0047 


0.45 


0.81 


4 


P'i^) 


EOV 


6-16 


0.6489 


0.0038 


1.98 


0.05 


4 


P%^) 


EOV 


8-16 


0.6449 


0.0031 


0.92 


0.47 


4 


P^i^) 


EOV 


7-16 


0.6498 


0.0029 


0.64 


0.70 


4 




EOV 


9-16 


0.6584 


0.0057 


1.28 


0.28 


4 




EOV 


6-16 


0.6488 


0.0036 


0.58 


0.77 


4 


P i^) 


EOV 


8-16 


0.6505 


0.0029 


0.48 


0.79 


4 




EOV 


8-16 


0.6568 


0.0041 


0.77 


0.57 


4 




EOV 


9-16 


0.6610 


0.0051 


1.24 


0.29 


4 


pi^) 


EOV 


7-16 


0.6524 


0.0035 


2.95 


0.07 


4 


A* 


point 


10-15 


0.926 


0.028 


0.17 


0.84 


4 


A T 

iV 


EOV 


2-15 


0.949 


0.010 


1.83 


0.10 


o 

8 


iV 


c 


8-15 


0.979 


0.008 


1.40 


0.23 


4 


iV' 


point 


NP 


NP 


NP 


NP 


NP 


4 


iV' 


EOV 


2-15 


1.289 


0.078 


1.83 


0.10 


8 


N' 


C 


8-15 


1.137 


0.089 


1.40 


0.23 


4 


A 


EOV 


6-15 


1.035 


0.010 


0.33 


0.92 


4 


A' 


EOV 


6-15 


1.302 


0.070 


0.33 


0.92 


4 



Table 3: Hadron masses aniq = 0.025. Notation: superscript is dimension of representation 
of time slice group; number of links in parenthesis; tilde (~) state has extra 70; notation 
abbreviated when unambiguous. 
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anig = 0.025 


particle 


source 


range 


mass 


error 


xVd-o.f. 


confidence 


parameters 


^* 

TT 


point 


7-16 


0.845 


0.072 


1.41 


0.21 


4 


^* 
TT 


EOV 


1-16 


0.842 


0.012 


1.13 


0.33 


4 


^* 
TT 


C 


1-16 


0.853 


0.006 


0.45 


0.94 


4 


/o/oo 


point 


9-16 


0.696 


0.007 


1.43 


0.22 


4 


/o/oo 


EOV 


8-16 


0.699 


0.010 


1.05 


0.39 


4 


/o/oo 


C 


8-16 


0.697 


0.015 


1.40 


0.23 


4 


ao(l) 


EOV 


8-16 


0.848 


0.030 


1.20 


0.30 


4 


«o(2) 


EOV 


6-16 


0.827 


0.018 


1.36 


0.22 


4 


ao(3) 


EOV 


6-16 


0.829 


0.020 


0.53 


0.82 


4 


ai 


point 


NP 


NP 


NP 


NP 


NP 


4 


ai 


EOV 


6-16 


0.886 


0.022 


1.38 


0.21 


4 


ai 


C 


6-16 


0.892 


0.021 


0.60 


0.76 


4 


a?(l) 


EOV 


7-16 


0.887 


0.036 


1.74 


0.11 


4 




EOV 


6-16 


0.905 


0.018 


0.93 


0.48 


4 


a?(2) 


EOV 


6-16 


1.023 


0.040 


0.58 


0.77 


4 


«?(2) 


EOV 


8-16 


0.927 


0.044 


0.77 


0.57 


4 


ai(3) 


EOV 


NP 


NP 


NP 


NP 


NP 


4 




point 


NP 


NP 


NP 


NP 


NP 


4 


w 


EOV 


7-16 


0.823 


0.067 


1.33 


0.24 


4 




C 


8-16 


1.011 


0.096 


1.00 


0.42 


4 




EOV 


6-16 


0.973 


0.044 


0.99 


0.43 


4 




EOV 


7-16 


0.843 


0.026 


0.95 


0.46 


4 




EOV 


7-16 


0.814 


0.050 


1.72 


0.11 


4 


m) 


EOV 


7-16 


0.855 


0.032 


0.65 


0.69 


4 


W{3) 


EOV 


7-16 


0.873 


0.063 


1.96 


0.07 


4 



Table 4: Hadron masses aniq = 0.025. Notation: superscript is dimension of representation 
of time slice group; number of links in parenthesis; tilde (~) state has extra 70; notation 
abbreviated when unambiguous. 
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/3 = 


5.85 










aniq = 


= 0.01 






aniq = 


0.025 




particle 


range 


mass 


error 


particle 


range 


mass 


error 


TT 


7-14 


0.2743 


0.0005 


TT 


9-14 


0.4243 


0.0008 


TT 


7-15 


0.4385 


0.0080 


TT 


6-14 


0.5577 


0.0048 


P 


6-14 


0.6476 


0.0149 


P 


6-13 


0.7183 


0.0056 


P 


4-12 


0.6258 


0.0084 


P 


4-12 


0.7126 


0.0040 


/o/oo 


5-13 


0.5624 


0.0321 


/o/oo 


6-14 


0.8075 


0.0314 


ai 


2-9 


0.8323 


0.0179 


ai 


4-12 


0.9832 


0.0284 


bi 


6-14 


1.553 


0.539 


bi 


6-13 


1.274 


0.213 


N 


5-13 


0.9501 


0.0276 


N 


5-13 


1.060 


0.008 


N' 


5-13 


0.7290 


0.0981 


N' 


5-13 


1.184 


0.065 



Table 5: Quenched liadron masses at /3 = 5.85. 









/3 = 


5.95 










aniq = 


= 0.01 






aniq = 


0.025 




particle 


range 


mass 


error 


particle 


range 


mass 


error 


TT 


6-14 


0.2501 


0.0009 


TT 


6-14 


0.3875 


0.0007 


TT 


8-16 


0.3215 


0.0044 


TT 


4-12 


0.4512 


0.0020 


P 


3-11 


0.5159 


0.0040 


P 


6-14 


0.5954 


0.0028 


P 


2-10 


0.5192 


0.0042 


P 


7-15 


0.5931 


0.0041 


/o/O 


8-16 


0.4777 


0.0541 


/o/O 


4-12 


0.6553 


0.0083 


ai 


2-10 


0.7184 


0.0090 


ai 


7-15 


0.8126 


0.0382 


bi 


3-11 


0.7073 


0.0228 


bi 


6-14 


0.8615 


0.0483 


N 


7-15 


0.7247 


0.0285 


N 


8-16 


0.8931 


0.0097 


N' 


7-15 


1.135 


0.220 


N' 


8-16 


0.9625 


0.115 



Table 6: Quenched hadron masses at /3 = 5.95. 





13 = 5.6 


aniq = 0.01 


/3 = 5.95 aniq = 0.0116 


particle 


mass 


error 


mass 


error 


TT 


0.2680 


0.0010 


0.2694 


0.0009 


TT 


0.351 


0.004 


0.3408 


0.0044 


P 


0.518 


0.004 


0.5244 


0.0040 


P 


0.515 


0.004 


0.5271 


0.0042 


N 


0.748 


0.004 


0.7426 


0.0285 



Table 7: Comparison between quenched hadron spectrum at /3 = 5.95 and aniq = 0.0116, 
and spectrum of full QCD at /3 = 5.6 and aniq = 0.01. 
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13 = 5.6 


anig = 0.025 


/3 = 5.95 anig = 0.029 


particle 


mass 


error 


mass 


error 


TT 


0.4189 


0.0005 


0.4173 


0.0008 


TT 


0.513 


0.006 


0.4810 


0.0020 


P 


0.637 


0.005 


0.6166 


0.0028 


P 


0.642 


0.004 


0.6128 


0.0041 


N 


0.981 


0.008 


0.9380 


0.0097 



Table 8: Comparison between quenched liadron spectrum at /3 = 5.95 and anig = 0.029, and 
spectrum of full QCD at /3 = 5.6 and anig = 0.025. 
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